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Abstract 

Based on Suzaku X-ray observations, we study the hot gas in regions around the NGC 4839 group of galaxies and 
the radio rehc in the outskirts of the Coma cluster From spectral analysis, the temperature of the gas shows a gradual 
decline from 5 keV around NGC4839 to about 3.6 keV at the radio relic. Across the radio relic, the temperature 
drops steeply by approximately a factor of 2 from 3.6 to 1 .5 ke V. This temperature drop can be interpreted as a shock 
with Mach number M = 2.2 + 0.5, which is also consistent with the surface brightness profile and corresponds to 
a shock velocity of Uj - (1410+ 110) km s We also derive a lower limit of B > 0.35// G on the magnetic field 
strength in the relic region from upper limits to inverse Compton X-ray emission from the non-thermal electron 
population that produce the radio relic. The existence of a shock front with magnetic fields suggests that it may be 
responsible for accelerating the non-thermal electrons. However, if they are accelerated according to the simplest 
theory of diffusive shock acceleration (DSA) in which test particles are injected from a thermal distribution, the 
electron spectrum expected from the measured Mach number and shock compression would be steeper than that 
inferred from the observed, spatially integrated spectrum of the radio relic with index a - 1.18, taking into account 
radiative loss effects. The expected electron acceleration efficiency may also be too low to be compatible with the 
observed radio flux. The radio relic may nevertheless result from DSA if it is initiated by the injection of a pre- 
existing population of non-thermal electrons. We confirm that the metal abundance around NGC 4839 is higher than 
in its vicinity and declines to a 0.2 solar value in the outermost regions. Combined with previous measurements, we 
discuss the temperature and abundance structure of the Coma cluster in the southwest direction. The temperature 
profile becomes shallower outside NGC4839 and is also slightly shallower in the outermost region. The metal 
abundance reaches a significant peak at NGC4839 and shows a radially decreasing trend toward the outskirts. Based 
on the correspondence between the temperature jump and the radio relic as well as the metal distribution, these facts 
can be interpreted in terms of a picture in which NGC 4839 is falling into the Coma cluster and the accompanying, 
metal-enriched gas is being stripped because of ram pressure. We suggest that the shock front associated with the 
radio relic is generated by the combined interaction of the pre-existing intracluster medium with gas conveyed by 
the infalling NGC 4839 group, together with the inflow of cooler gas from the large-scale structure filament that 
connects Coma with Abell 1367 in the southwest as part of the Coma supercluster 

Keywords: galaxies: clusters: individual (Coma, NGC4839) — galaxies: intergalactic medium — Shockwaves 
— X-rays: galaxies: clusters 



1. Introduction 

According to modern theories of large scale structure forma- 
tion, clusters of galaxies are built up hierarchically through suc- 
cessive mergers and accretion (Voit 2005; Kravtsov & Borgani 
2012). In particular, major mergers involving multiple subclus- 
ters of comparable masses (hereafter simply "mergers" unless 
otherwise indicated) are dramatic events that strongly affect the 
intracluster medium (ICM) and can convert as much as lO*"^ 
erg of kinetic energy into thermal energy through shocks and 
turbulence (Markevitch & Vikhlinin 2007). Shocks not only 

* Last update: 2013 February 13 



heat the gas but can also produce high-energy, nonthermal par- 
ticles through the mechanism of diffusive shock acceleration 
(DSA) (Blandford & Eichler 1987) and provide a key observa- 
tional tool in the study of cluster growth. X-ray observations 
have revealed direct evidence of shock fronts through discon- 
tinuities in the surface brightness and temperature in the ICM 
of merging clusters, for example the so-called "Bullet" clus- 
ter IE 0657-56 (Markevich et al. 2002), and more recently in 
AbeU 2146 (Russefl et al. 2010) and Abell 754 (Macario et al. 
201 1). Such shocks discovered so far in X-rays are associated 
with clusters undergoing relatively major mergers (Markevitch 
2010). Large-scale shocks driven by the accretion of diffuse 
gas that is not bound in viriahzed systems are also theoretically 
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Fig. 1. Left: Background subtracted ROSAT All Sky Survey image in the 0.1-2.4 keV band of the southwestern regions of the Coma cluster. The color 
scale unit is cts Ms^'pixel"^ (1 pixel=45" x45"). The white boxes show the FOVs of the Suzaku XIS observations we discuss, while the magenta circle 
shows Coma's virial radius r2oo. corresponding to 83.'9 from the cluster center Right: NXB subtracted Suzaku FI+BI mosaic image in 0.5-8.0 keV band, 
smoothed with a 2-dimensional gaussian with o" = 16 pixel =17" . The image is corrected for exposure time but not for vignetting. The color scale unit 
is cts Ms^'pixel"^ (1 pixel=l."04 X l."04). The white and gray dashed annuli show the regions used for the spectral analysis. Small green circles show 
point sources. The NVSS 1.4 GHz image of the radio relic is shown in yellow contours. 



Table 1. Observation log and exposure time after data screening without and with COR-cut (COR2 > 6 GV) 



Name (Obs I.D) 


(R.A., DEC) 


Observation Date 


Exp. Xime' 


Exp. Xime* 


45' (802047010) 


(194.26, 27.58) 


2007-12-02X10:08:56 


28.0 


22.5 


60' (802048010) 


(194.03,27.44) 


2007-12-04X04:47:49 


32.9 


27.6 


Radio relic (803051010) 


(193.86, 27.30) 


2008-12-23X21:06:23 


172.3 


145.9 


SW55 (806048010) 


(193.64, 27.07) 


2011-06-23X20:47:37 


15.2 


14.6 


SW60 (806049010) 


(193.53,26.98) 


2011-06-24X07:17:49 


12.7 


8.8 


COMABKG (802083010) 


(198.75,31.63) 


2007-06-21X00:14:27 


28.2 


24.1 


t: COR2 > GV 


COR2 > 6 GV 









expected and generally seen in numerical simulations (Miniati 
et al. 2000; Ryu et al. 2003), but have yet to be observationally 
confirmed in X-rays. 

Xhe Coma cluster is one of the best-studied clusters of galax- 
ies because of its low redshift (z = 0.023) and its high richness 
(Biviano 1998). In view of its inhomogeneous distribution of 
galaxies and disturbed X-ray surface brightness, it is acknowl- 
edged to be a merging system (Hughes et al. 1988; Hughes et 
al. 1993). Xhe ROSAX X-ray image clearly shows some sub- 
structures (Briel et al. 1992). Watanabe et al. (1999) presented 
temperature maps based on ASCA data that revealed the ICM 
distribution to be non-isothermal. Neumann et al. (2001) re- 
ported the structure of Coma on large scales based on XMM- 
Newton mosaicking observations, which indicate that it is now 
undergoing at least one subcluster-cluster merger In particu- 
Im; the group of galaxies associated with the elUptical galaxy 



NGC 4839 is a prominent substructure located southwest of 
the cluster center, and may be falling into Coma (Adami et al. 
2005). Neumann et al. (2001) also found indications of a bow 
shock and ram pressure stripping around NGC 4839. 

Merging clusters including Coma also often exhibit radio 
halos and/or radio relics, which are spatially-extended, syn- 
chrotron emission from non-thermal populations of electrons 
that are somehow accelerated over scales of up to ~Mpc 
(Ferrari et al. 2008; Feretti et al. 2012). Radio halos oc- 
cur in the central regions with a fairly symmetric morphol- 
ogy, and are often considered to result from stochastic ac- 
celeration by merger-induced turbulence (e.g., Brunetti 201 1). 
Contrastingly, radio relics are found near the cluster periph- 
ery with a tangentially-elongated morphology, and may be 
due to DSA in accretion or merger shocks (Ensslin et al. 
1998;Bruggen et al. 2012). 
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The Coma radio relic is located near its virial radius about 2 
Mpc from the cluster center, beyond NGC4839 but in the same 
southwest direction. Feretti & Neumann (2006) investigated 
the relic region with XMM-Newton observations. Although 
they detected thermal emission around the radio relic, they 
did not find any temperature enhancement or evidence of a 
shock in regions northeast of the relic towards the cluster cen- 
ter. Thus they proposed that the relativistic electrons are accel- 
erated by turbulence generated by the infall of the NGC 4839 
group, rather than by shock acceleration. Brown & Rudnick 
(2011) discussed the possible connection with ongoing accre- 
tion based on their radio and optical observations. Recently, 
Planck found signatures of shock fronts through measurements 
of the Sunyaev-Zeldovich effect in two regions about half a 
degree to the west and to the south-east of the cluster cen- 
ter, which corresponds to the edge of the radio halo (Planck 
Collaboration et al. 2012). They derive Mach numbers Alw = 
1.95!°;;^ and Mse = 2.03+°-,'4 for the shocks in the west and 
southeast, respectively. 

Possible associations between shocks and radio relics have 
also been found in a number of other merging clusters, van 
Weeren et al. (2010) reported the discovery of a double radio 
relic with an unprecedentedly narrow morphology in a merging 
cluster. The spatially resolved radio spectral index is c ~ -0.6 
along the sharp outer edges of the relics, which presumably re- 
flects the spectrum of the emitting electrons that have just been 
accelerated and are yet to be afffected by radiative losses. Under 
the assumption that these electrons are generated according to 
the simplest, test-particle theory of DSA, the spectral index 
can be related to the Mach number of the shock, which in this 
case would be Al ~ 4.5. Subsequent, deep X-ray follow-up ob- 
servations revealed the existence of a shock front at the radio 
relic (Ogrean et al. 2012; Akamatsu & Kawahara 2011). The 
cluster Abell 3667 has a pair of radio relics, of which the north- 
western one is the brightest diffuse cluster radio source in the 
sky (Rottgering et al. 1997). Finoguenov et al. (2010) found 
evidence of a merger-induced shock at the outer edge of this 
relic, and showed that it was energetically capable of powering 
the radio relic. Thus, it appears that radio relics can be good 
tracers of merger shocks in clusters. 

In this paper, we discuss X-ray observations focusing on 
the regions of interaction between Coma and NGC 4839 and 
the radio relic located near the outskirts of the Coma cluster. 
We assume cosmological parameters Hq = 70 km s"' Mpc"', 
Q.M = 0.27 and Qa - 0.73, which implies 28.9 kpc per ar- 
cminute at Coma's redshift of z = 0.023. The virial radius is 
approximately raoo = 2.77h-l^({T)/10keWy'^/E(z) Mpc, where 
E(z) = (Qm(1 +z)^ + 1 - ^m)^'^ (Henry et al. 2009). For our 
cosmology and average temperature k{T} = 7.8 keV as used 
in Matsushita (2011), r2oo = 2.43 Mpc, corresponding to 83.'9 
from the cluster center As our fiducial reference for the solar 
photospheric abundances, we adopt Lodders (2003). Unless 
otherwise stated, the errors correspond to 68% confidence for 
a single parameter 

2. Observations and Data Reduction 

As shown in Fig. 1, Suzaku carried out five observations with 
different pointings of the southwest region of the Coma clus- 



ter in December 2007, 2008 and July 201 1. They covered the 
regions from around NGC4839 to beyond the radio relic. We 
refer to these pointings as 45', 60', Radio Relic, SW55 and 
S W60. The observation log is given in Table 1 . All the obser- 
vations were performed with either normal 5 x 5 or 3 x 3 clock- 
ing modes. The virial radius r2oo of the Coma cluster (2.43 
Mpc corresponding to 83. '9) is indicated by the magenta cir- 
cles in Fig. 1. 




0.5 1 2 

Energy (keV) 



Fig. 2. The spectram of the COMA BGK field used for the background 
estimation, after subtraction of the NXB and point sources. The XIS BI 
(black) and FI (red) spectra are fitted with CXB + Galactic components 
(LHB, MWH) (apec+wabs(apec+powerlaw)). The CXB spectrum is 
shown with a cyan curve, and the LHB and MWH components are 
indicated by green and blue curves, respectively. 

The XIS instrument consists of 4 CCD chips: one back- 
illuminated (BI: XISl) and fliree front-illuminated (FI: XISO, 
XIS2, XIS3) detectors. We used HEAsoft version 6.12 and 
CALDB 2012-02-10 for all of the Suzaku data analysis pre- 
sented here. We performed event screening with the cosmic- 
ray cut-off rigidity COR2 > 6 GV to increase the signal to noise 
ratio. We also limited the Earth rim ELEVATION to > 10° to 
avoid contamination from scattered solar X-rays from the day- 
side Earth limb. For the SW55 and SW60 data, we applied 
additional processing for the XIS 1 detector to reduce the NXB 
level, which increased after the change in the amount of charge 
injection in 2010 June. The detailed processing procedures are 
the same as those described in XIS analysis topics^. 

We extracted pulse-height spectra in 12 annular regions 
whose boundaries are given in Table 3 with the center at 
(12''57"24\27°29'54"), which is the X-ray peak associated 
with NGC4839 (Neumann et al. 2001). For the 43.'0- 54.'0 
annulus, we combined the SW55 and SW60 data to increase 
the signal to noise ratio. There is a missing annulus from 
38.'0- 43.'0. This region was only covered by the SW55 ob- 
servation, and the shorter exposure and lower surface bright- 
ness led to too few counts for spectral fitting. In all annuli, the 
positions of the calibration sources were masked out using the 
calmask calibration database (CALDB) file. 

' http://www.astro.isas.jaxa.jp/suzaku/analysis/xis/xisl_ci_6_nxb/ 
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NOMINAL 


CXBMIN-3%ofNXB 


CXBMAX + 3% of NXB 


CONTAM-10% 


CONTAM-i-10% 


LHB 












kT (keV) 


0.08 (fix) 


0.08 (fix) 


0.08 (fix) 


0.08 (fix) 


0.08 (fix) 


norm* (xlO"^) 


9.7 ±3.2 


11.8 


7.2 


6.9 


14.9 


MWH 












kT (keV) 


0.18 ±0.02 


0.19 


0.14 


0.15 


0.19 


norm* (xlO""*) 


9.8 ±2.7 


9.7 


17.7 


15.1 


9.9 


CXB (r= 1.41 fix) 












norm' (xlO""*) 


8.7+0.4 


7.8 (fix) 


9.5 (lix) 


8.7 (lix) 


8.7 (lix) 




63/82 


66/83 


73/83 


63/83 


63/83 



*: Normalization of the apec component scaled by a factor l/4007r. 

Nonn=^gL_ J^^^^jy/(4;j-(l +z^)D^) X lO"^'' cm"^ arcmin"-^, where Z)a is the angular diameter distance to the source, 
t: The CXB intensity normalization in Kushino et al. 2002 is 9.6 X 10""* for T = 1.41 
in units of photons keV"' cm"^ s ' at 1 keV. 



Table 3. Best-fit spectral parameters of the ICM 



Region* 



0' 


-3' 


3.'0 


-5.'5 


5. '5 


-10' 


10' 


-13' 


13' 


-16' 


16' 


-19' 


19' 


-22' 


22' 


-26' 


26' 


-30' 


30' 


-34' 


34' 


-38' 


43' 


-54' 



kr (keV) 

-'•■-'^-0.26 

4 85+"-^^ 

;+0.20 
-0.18 

3^74+0.21 



Z (Zq) norm' /d.o.f 



0.60 
0.63 



+0.10 
-0.09 
+0.13 



3.76+ 



^-0.16 
+0.17 
0.16 

+0.27 



3.64 

3-47-0.25 
3.64!«:g 

T 40+0.16 
^•^^-0.14 

T ^7+0.21 
-'-"-0.20 

^•"'-0.22 
1 54+0.45 



-0.13 

Q 44+0.12 

^•^^-0.13 
^•-"-0.09 

0-36_o.o8 
n ^/;+0 09 
0-36_o.o8 

o.29!;j:|^ 

0.28-- 

2r°°^ 

^■^ '-0.07 

17+°- 12 

^•^ -0.09 

0.15 (fix) 



+rr 

-2.0 
+2.0 

43.9^? 

i:8 

-1.5 
+ 1.3 
-1.0 
+0.9 
1.0 
+ 1.3 
0.7 
+0.4 
-0.3 
+0.3 
0.3 
+0.3 
-0.3 
+0.6 
0.4 
+0.5 
0.3 



60.9 
50.4 



39.1 

29.7 

23.7 

17.2 

14.2 

9.4 

7.9 

5.5 

2.3 



183/154 
111/129 
105 / 125 
143/137 
185/170 
183/175 
114/110 
402/360 
395 / 360 
361/360 
208/ 174 
118/147 





North Direction (Gray dotted annuli in Fi; 


g- 1) 




3.'0 

5. '5 


-5'5 7 77+i)-'9 1 5+0"' 30 4+'' 1 

'•''-0.71 "-'^-O.H -'"-^-1.5 ' 

-10' 6 50+"-^'^ 30^0-'^ 75 7+"'' 1 

i\i u.JVJ_,)^45 u.JVJ_,)^j4 1 


28/ : 

20 / . 


112 
112 



*: Radius from NGC4839, 

t:Normalization of the apec component scaled with a factor l/400;r. 
nona=^ ^ nenHdVI(An(l +:^)D\) x IQ-^" cm-^ arcmin-^, 
where Da is the angular diameter distance to the source. 



3. Spectral Analysis and Results 

3.1. Method of Spectral Fits 

The observed spectrum was assumed to consist of optically 
thin thermal plasma emission from the ICM, from the Local 
Hot Bubble (LHB) and the Milky Way Halo (MWH) as the 
Galactic foreground components, cosmic X-ray background 
(CXB), and non-X-ray background (NXB). The NXB com- 
ponent was estimated from the dark Earth database using the 
xisnxbgen program in FTOOLS (Tawa et al. 2008), and was 
subtracted from the data before the spectral fit. To adjust for 
the long-term variation of the XIS background due to radiation 
damage, we accumulated the NXB data for the period between 
150 days before until 150 days after each observation. 



For the spectral fits, we used XSPEC verl2.7.0. Using the 
xissimarfgen program in FTOOLS (Ishisaki et al. 2007), we 
generated spatially uniform ARFs over a circular region of 
20.'0 radius for the background and ICM emissions. We car- 
ried out spectral fits to the pulse-height spectrum in each an- 
nulus separately. We analyzed the spectra in the 0.35-7 keV 
range for the BI detectors and 0.5-8 keV for the FT detectors. 
Because of possible contamination by solar wind charge ex- 
change, we only analyzed the spectra above 0.9 keV for all de- 
tectors in the SW55 region (43 .'0- 54.'0 annulus). None of the 
other spectra were significantly afi'ected by solar wind charge 
exchange. In addition, above 5 keV, the CXB dominates the 
observed spectra of the 43.'0- 54.'0 annulus. For this region, 
we only use the spectrum below 5 keV to increase the signal 
to noise ratio. As shown by the green circles in Figure 1, we 
found 6 point-like sources and masked out the pixels within a 
1 arcmin radius around these sources. 

In order to estimate the CXB and Galactic background com- 
ponents, we selected the ofi'set Suzaku observation of COMA 
BGK (Observation ID = 802083010), which was located at 
6 degrees angular distance from NGC 4839. Although the 
background components were estimated in Takei et al. (2008) 
and Sato et al. (2011), we re-analyzed the same data to de- 
termine the systematic errors. We modeled the sky back- 
grounds as a combination of the Local Hot Bubble (LHB ~ 
0.08 keV), the Milky Way Halo (MWH ~ 0.3 keV), and the 
cosmic X-ray background (CXB), using the spectral model: 
apec + wabs{apec + powerlaw). In the fits, we fixed the tem- 
perature of the LHB component to 0.08 keV. The redshift and 
metal abundance of both the apec components were fixed at 
and one solar, respectively. The resultant spectrum and pa- 
rameters are shown in Figure 2 and Table 2. The temperature 
of the MWH is 0.18 ± 0.02 keV. The temperature and intensity 
are consistent with the typical Galactic emission and previous 
measurements (Yoshino et al. 2009; Takei et al. 2008; Sato et 
al. 2012). For the estimate of the systematic error of the back- 
ground spectrum, we adopted the ampUtude of the CXB fluctu- 
ations as 10%, the uncertainty of the NXB intensity as 3% and 
the uncertainty in the contamination on the optical blocking 
filter to be ±10%, respectively. Using this background model, 
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(a) 0'-3' 



(b)3'-5.'5 



(c)5.'5-10' 





(d) 10'-13' 



Energy (keV) 



(e) 13'-16' 



Energy (keV) 



(f) 16'-19' 



Energy (keV) 




(g) 19'-22' 



Energy (keV) 



(h) 22'-26' 



Energy (keV) 



(i) 26'-30' 



1 2 
Energy (keV) 




Energy (keV) 



Energy (keV) 



Energy (keV) 



Fig. 3. NXB subtracted spectra in eacli annular region. The XIS BI (Black) and FI (Red) spectra are fitted with the ICM model (Magenta: wabs X apec), 
along with the sum of the CXB and the Galactic emission {apec + wabs(apec + powerlaw)). The CXB component is shown with a black curve, and the 
LHB and MWH emissions are indicated by green and blue curves, respectively. For the 43'-54' annulus, we show the SW60 spectrum. 



we also examined the effect of systematic errors on our spectral 
parameters (Sec. 3.3). 

The two normalizations and temperatures of the MWH com- 
ponent, and the normalization of the power-law model for the 
CXB component were allowed to vary within the range of the 
errors in the background estimate (Table 2). We employed a 
single temperature thermal model (wabs x apec) for the ICM 



emission of the Coma cluster In the radio relic region, in addi- 
tion above model, we employed two thermal model to taken ac- 
count for an effect of PSF mixing and projection (see sec. 3.2). 
The interstellar absorption was fixed using the 21 cm measure- 
ment of the hydrogen column, A^h = 1.0 x 10^° cm"^ (Dickey 
& Lockman 1990). The redshift of the ICM component was 
fixed at 0.023 (Biviano 1998). In the inner region (r < 38' 
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from NGC4839), the temperature, abundance, and normaliza- 
tion of the single-temperature model were free parameters. In 
the outer regions (r > 38'), we fixed the ICM metal abundance 
to 0.15. In the simultaneous fit of the BI and FT data, only 
the normalizations were allowed to differ, although we found 
that the derived normahzations were consistent within 15%. In 
general, each data set was fitted well with the above mentioned 
model. 

3.2. Results 

Temperature and abundance are the basic parameters that we 
derive from the spectral analysis of the ICM. In order to inves- 
tigate the radial trend of the temperature and metal abundance, 
we analyzed the X-ray spectra in annuli as described in Sec. 2 
and Galactic absorption with A^h = 1-0 x 10^" cm"^ (Dickey 
& Lockman 1990). . The spectra and the best-fit models 
for all the annuh are shown in Fig. 4. The parameters and 
the resultant values are listed in Table 3. The temperature 
and abundance within 3' from NGC4835 are kT = 4.18 + 0.17 
keV and Z - 0.60 ±0.10 Zg, respectively, which is are con- 
sistent with the XMM-Newton results of kT = 4.4+^-^ keV 
and Z = 0.36^q }q Z© (Neumann et al. 2001). The tempera- 
ture profile is almost flat at kT ~ 3.6 keV between r = 10' and 
r = 30', which corresponds to 500 kpc. Those values are also 
in good agreement with previous measurements (Neumann et 
al. 2001; Sato etal. 2011). 

The abundance profile shows a peak at NGC4839 and grad- 
ually decreases with radius from 0.6 to 0.17 solar. The high 
abundance value around NGC4839 is also consistent with the 
previous XMM-Newton results. In contrast, the ICM to the 
north of NGC4839 (gray dotted annuli in Fig. 1) has differ- 
ent properties compared with same radii to the southwest. The 
ICM temperature is higher than beyond NGC4839, as has al- 
ready been reported by Neumann et al. (2001). They inter- 
preted this temperature structure as a bow shock associated 
with the infaU of NGC4839. The abundance is smaller than 
to the SW, which may indicate that the SW gas is affected by 
metal enrichment from NGC4839 via ram-pressure stripping. 
We will discuss the enrichment process in the cluster outskirts 
in Sec. 4.2. In addition, there is a small deviation from previous 
work (Sato et al. 201 1). The deviation is possibly explained by 
the difference in the extraction areas used for the spectral anal- 
ysis. 

Previous X-ray observations did not find evidence of a shock 
associated with radio relic in the Coma cluster (Neumann et al. 
2001; Feretti & Neumann 2006). The jump in our temperature 
profile suggests that there is a shock front in the relic region. 
The properties of the shock front will be discussed further in 
Sec. 4.3. 

The shock front expected to be located at the edge of the 
reUc. However, in our analysis, the high temperature com- 
ponent dose not show the inner edge of the relic. Previous 
our studies of radio relic clusters (A3667, A3376: Akamatsu 
et al. 2012; Akamatsu et al. 2012) showed a possibility multi- 
temperature plasma in those region. It is likely that this came 
from a mixing effect due to large PSF of Suzaku XRT and pos- 
sibly projection effect. To examine the those effect in the radio 
relic region (34'-37'), we included additional thermal compo- 
nent (2kT) to the single temperature model. We obtained an al- 



most acceptable fit withx^/ d.o.f . = 204/172 compared with the 
single temperature case of 208/174. The resultant hot compo- 
nent temperature was derived as kThigi, = S.Sb^Q^'^ keV, and the 
cool one as kTiow=^ -^O^qIq keV while the common abundance 
was fixed to 0. 15 solar. Those values are almost consistent with 
single temperature fits of inside (30' - 34' : kT = 3.57;|;g 2Q) and 
outside (43' -54':kT = l .54+° ^5) of the radio relic. 

3.3. Systematic Errors in the Radial Spectral Analysis 

In order to discuss the ICM properties in the cluster outskirts, 
systematic errors must be carefully addressed. We considered 
three components for the systematic error: the NXB estimation 
uncertainty, fluctuations in the CXB intensity, and uncertain- 
ties in the contamination on the optical blocking filter. The 
NXB uncertainty was estimated by Tawa et al. (2008) as +3%. 
To estimate the fluctuation of the CXB, we follow the method 
used in Hoshino et al. (2010). They scaled the fluctuation mea- 
sured by Ginga (Hayashida et al. 1989) with the field of view 
and the flux detection limit. We removed point sources identi- 
fied by eye in the Suzaku images as shown by the green circles 
in Fig. 1. To estimate fluctuations in each annular region, we 
adopt a conservative flux limit, 5e = 1 x lO"'-' erg s"' cm"^. 
Table 4 shows the resultant CXB fluctuations. Thus, the esti- 
mated CXB fluctuations are the range of 15-30 %. 

Another potentiall large systematic error component comes 
from contamination on the blocking filter of the XIS instru- 
ment. After the launch of Suzaku in 2005, the contamination 
kept increasing for two years, after which it became stable. 
^. All observations analyzed here were performed after 2007. 
The amount of contamination is now known with fairly high 
accuracy, for which we assign a 10% error. In addition to the 
above systematic errors, we also examined the consequences of 
adopting a different reference for the solar abundance, namely 
Anders & Grevesse (1989) instead of Lodders (2003). 

The resultant range of best fit parameters after accounting 
for the systematic errors are shown in Fig. 4 as dotted lines 
in blue for CXBh-NXB, green for filter contamination and ma- 
genta for different solar abundances. The effects of systematic 
errors become larger toward the outskirts. For most regions , 
the systematic errors are smaller or comparable to the statis- 
tical errors. In the outermost region, the resulting systematic 
uncertainty in the ICM temperature is about 0.5 keV. In the 
abundance profile, the difference between the two abundance 
tables is significant around NGC4839. The difference in the 
abundances fit with the Lodders (2003) and Anders & Grevesse 
(1989) models is close to the difference in the iron abundance 
between the two tables. Thus, the fitted abundances are very 
similar in absolute terms 

4. Discussion 

Suzaku performed five observations with different pointings 
of the southwestern regions of the Coma cluster. The radial 
profile of the temperature was measured beyond the radio reUc 
region as well as beyond r2oo (~ 83.'9 from the center) for 
the first time. The ICM temperature declines from ~5.0 keV 
around NGC4839 to ~3.6 keV at the radio relic. Furthermore, 

^ http://www.astro.isas.jaxa.jp/suzaku/doc/suzaku_tcl/nodel0.html#contami 
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Fig. 4. Radial profiles of the (a) ICM temperature, and (b) metal abundance, centered on NGC4839. The Suzaku best-fit values with statical errors are 
shown with red crosses. Gray crosses show the results for the northern regions (the gray dotted annuli in Figure 1). The gray dashed vertical lines show 
the approximate radial boundaries of the radio relic (1253+275). The black crosses show the previous measurements by Sato et al. (201 1), with converted 
1(T errors. The range of unceitainties due to the combined 3% variation of the NXB level and the maximum/minimum fluctuation in the CXB is shown 
by two blue dashed curves. In a similar way, the range of uncertainties due to the contamination on the blocking filter is shown by green dashed curves. 
Magenta dashed lines show changes due to different abundance tables. 



Table 4. Estimation of the CXB fluctuation. 



Region* 0'-3' 3.'0-5.'5 5.'5-10' 10'-13' 


13'-16' 


16'-19' 19'-22' 22'-26' 


26'-30' 


30'-34' 


34'-38' 


43'-54' 


Area' 28.3 60.9 95.9 55.7 
Flue* 37.0 25.2 20.1 26.4 


73.0 
23.0 

r 


88.7 79.1 94.8 
20.9 22.1 20.2 


125.0 
17.6 


115.9 

18.3 

-ii - 


109.9 
18.8 

1 -2 rr» 1 


145.1 
14.2 



*: Radius from NGC4839, f: Observational area [arcmin'^], p. CXB fluctuation assuming 5^. = 1x10 erg s ' cm ^ [%]. 



we find a steep drop in temperature from 3.6 to 1.5 keV be- 
yond the radio relic, which suggests the presence of a shock 
front. The metal abundance also shows a high value of Z ~ 0.6 
solar at NGC4839 and gradually decreases to 0.15 solar around 
the relic. Below, we discuss the ICM properties and their im- 
plications. 

4.1. Temperature Structure 

The radial profile of the ICM temperature provides infor- 
mation about the thermal and merger history of the cluster 
Previous studies of large cluster samples with ASCA and 
BeppoSAX have indicated that temperature decUnes with ra- 
dius out to at least ~ 0.4 raoo (Markevitch et al. 1998; De 
Grandi & Molendi 2002). Those results were extended to 
0.7 ^200 by more recent observations with Chandra and XMM- 
Newton (Vikhlinin et al. 2005; Pratt et al. 2007). Suzaku can 
extend temperature profiles up to rioo because of its low and 
stable instrumental background (George et al. 2008; Reiprich 
et al. 2009; Bautz et al. 2009; Hoshino et al. 2010; Simionescu 
et al. 2011; Sato et al. 2012; Walker et al. 2012). Previous re- 
sults indicate that the temperature in relaxed clusters shows a 



systematic drop by a factor of 2-4 from near the center to raoo. 

Here, we discuss the temperature structure in Coma, 
including previous XMM-Newton (Matsushita 2011) and 
Suzaku (Sato et al. 201 1) measurements. The profile is shown 
in Figure 5(a). The temperature was measured out to r2oo 
(83.'9), and shows a drop by a factor of 4 from the center to 
'"200- The amplitude of this drop is consistent with other clus- 
ters. To quantify the temperature drop, we fit these temperature 
profiles with a simple linear model. 



T{r)^{T)\a + b 



'"200 



(1) 



where we take k{T) - 7.8 keV as discussed in Sec. 1. We also 
consider the polytropic function 



T{r) = To 



1+1^ 

rr 



-1.5/3(7-1) 



(2) 



where Tq is the fitted central temperature of the cluster 
The polytropic model was used to fit a large cluster sam- 
ple in Markevitch et al. (1998). Previous studies with 
ASCA (Markevitch et al. 1998) report y = 1 .24+[j-J^ (90% error) 
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Fig. 5. Scaled radial profiles of (a) the noiTnalized ICM temperature, and (b) the metal abundance. The Suzaku best fit values with statical errors are 
shown with red crosses. Black and gray crosses show the data from the Virgo and Perseus clusters (Urban et al. 2011; Simionescu et al. 2011), and the 
blue crosses indicate the Coma results from XMM-Newton (Matsushita 2011). To normalize temperature profiles, we adopted average temperatures of 
7.8, 2.3 and 6.4 keV for the Coma, Virgo and Perseus clusters, respectively. We estimate ('200 from the average temperature (see Sec. 1). In (a), the blue 
and red dotted fine show the best fit Hnear functions for the XMM-Newton and Suzaku data, respectively. (See the text for the details.) In (b), all of the 
abundances are converted to the Lodders (2003) solar values. 



Table S. Best-fit parameters for the temperature profiles 





Linear 


Polytrophic 




a 


b 


y 


Markevitch et al. (1998) 








Pratt et al. (2007) 


1.19 


-0.74 




This work (0.0 - 0.5 noo) 


0.99 ± 0.04 


-0.76 ± 0.03 


1.25 ±0.04* 


This work (0.5 - 0.9 r^oo) 


0.72 ±0.06 


-0.39 ±0.08 


1.27 ±0.03"^ 



*:Fitting range is 0.1 - 0.5 r^oo- 



t: 7*0 set to same value determined in inner region. 

in 30 clusters. We adopt typical values for the gas density pro- 
file of p - 0.6 and r^. - 9' (Cavaliere & Fusco-Femiano 1978). 
We fit the temperature profile in the inner (0.1 - 0.5 r2oo) and 
outer (0.5 - 0.9 r2oo) regions with these two models. The re- 
sulting values are summarized in Table 5. 

There is a clear difference in a temperature gradient between 
the inner and outer regions. For inner region, our best-fit val- 
ues of the gradient b = -0.76 + 0.03 and parameter y=1.25 + 
0.03 is almost agree with previous studies (Markevitch et al. 
1998; Finoguenov et al. 2001; Pratt et al. 2007). On the other 
hand, in the outer region the temperature gradient b is shal- 
lower (b = -0.39 + 0.08). Because the polytrophic model is 
sensitive to the inner structure, it is not surprising that the y pa- 
rameter does not show any diff'erence. The flatter temperature 
profile beyond 0.5 r2oo indicates that the outer gas has under- 
gone some additional heating. The effect of heating related to 
merger activity on the ICM temperature is confirmed in many 
systems (Markevitch & Vikhlinin 2007). Although we studied 
a very limited area of the Coma cluster, it is likely that south- 
west of Coma has experienced heating associated with the in- 
fall of the NGC4839 group. In order to explore the effects of 



more continuous accretion, we need temperature profiles in di- 
rections which are not affected by recent merger activity. 

4.2. Abundance Structure 

Measurements of the abundances in the outskirts of clus- 
ters are important to study the origin of metals in the ICM. 
Recently, metal abundances at large radii were reported in 
A399 and A401, and in the Perseus, Virgo and Hydra clus- 
ters (Fujita et al. 2008; Simionescu et al. 2011; Urban et al. 
2011; Sato et al. 2012). These abundance profiles extend 
out to r2oo, where the metallicity was found to be 0.1-0.3 so- 
lar Previously, several metal enrichment mechanisms have 
been proposed (e.g., Schindler & Diaferio 2008). One is ram- 
pressure stripping; the metal enriched gas in the galaxies is 
stripped by ram pressure when the galaxies move through the 
ICM (Gunn & Gott 1972; Fujita & Nagashima 1999). HI 
observations have provided many examples of galaxies being 
stripped by ram-pressure. In particular, the Virgo cluster has 
several spiral galaxies affected by ram-pressure stripping (e.g. 
NGC4522: Kenney et al. 2004; NGC4402: Crowl et al. 2005; 
NGC4472: Kraft et al. 201 1). X-ray observations by Sun et al. 
(2010) found hot tails associated with infalling galaxies onto 
A3627 (ESO 137-001 and ESO 137-002). Those stripping 
events are important not only for galaxy evolution but also for 
the abundance evolution of the ICM. 

Matsushita (2011) reported radial profiles of the Fe abun- 
dance out to 0.3-0.5 r2()o for 28 clusters of galaxies observed 
with XMM-Newton. These clusters show similar abundances 
within 0.3 r2oo of about 0.4-0.5 solar, and a decrease to -0.3 
solar at 0.5 r2oo, using the Lodders (2003) abundance table. 

In the Coma cluster, our abundance profile shows a sig- 
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nificant peak at 0.5r2oo which corresponds to the position 
of NGC4839. Also, abundances outside of NGC4839 (r > 
0.5 r2oo) are higher than or comparable to the central region 
(r < 0.2 r2oo)- At cluster outskirts which are fart from the clus- 
ter, the contributions from the central region (e.g., cD galaxy) 
are small. In addition, abundances to the SW of NGC4839 
(red crosses in Fig. 4) are higher than those to the north (grey 
crosses in Fig. 4). Since it is likely that NGC4839 in falling 
relative to the ICM from the SW to the NE, and there is also 
an increased density of gas to the SW of this galaxy, it is likely 
that the gas to the SW has been stripped from this galaxy or 
its associated group. Our data provide direct observational evi- 
dence that the ICM is enriched by infalling galaxies and groups 
at large distances from cluster center (0.5-0.8 rioo). 

Aside from the high abundance region associated with 
NGC4839, our data gives the abundance up to 0.9 r2oo. The ob- 
served abundance at large radii decreases to 0. 17 + 0. 1 1 solar. 
The value is similar to previously obtained abundance values 
at large radii in other clusters as mentioned above (Perseus and 
A399 & A401). Recently, a low abundance value was found 
at large radii in the Virgo cluster (Urban et al. 201 1). They in- 
terpreted the low abundance as arising from biasing by the Fe 
L Unes used to determine the abundance. The Virgo cluster is 
a low temperature system (k{T) ~ 2.5 keV) and the tempera- 
ture at r2oo has decreased to ~1 keV. At those temperatures, the 
metallicity would be determined primarily from the Fe-L Unes, 
which are sensitive to the temperature structure (Buote 2000). 
To examine the effect of the Fe L complex in our Coma spectral 
fits, we fitted the 30'-34' and 34'-38' annuU in the energy range 
1.5-8.0 keV for all the detectors. The temperature and normal- 
ization for each region were allowed to vary, but the abundance 
was required to be the same in both regions. The resultant ICM 
abundance was Z = 0.23 ± 0.09 Zq, which is consistent with 
our previous value (Z = 0.17 + 0.1 1 Zq). For the 1.5-8.0 keV 
energy range, the abundance would be determined by the Fe K 
lines. Thus, the abundance around 0.9r2oo in the Coma clus- 
ter is likely around 0.2 solar, and is not strongly affected by 
biasing from Fe L lines. 

We confirm the existence of metals in the cluster periphery 
in Coma along the direction of the merger with the NGC4839 
group. What does this tell us about the origin of these heavy 
elements? First, the metals might be due to ram pressure strip- 
ping from the NGC4839 group. The pattern of gas densities 
and abundances is consistent with this scenario. One concern 
is that the ICM densities at these large radii are small («e = 
6.1^Q^ X lO^^cm"-' around r2oo in our measurements), which 
might make ram pressure ineffective (Fujita & Nagashima 
1999). After the metals are shipped from tiie NGC4839 group, 
they could diffuse due to the turbulence associated with the 
merger. The outer regions of clusters might also be polluted 
with metals due to galactic superwinds, which may have op- 
erated at z ~ 3 (Ferrara et al. 2005) before cluster formation. 
Some numerical simulations (Cen & Ostriker 2006) suggest 
that galactic superwind feedback at high redshifts driven by 
active galactic nuclei and star bursts could provide metals to 
the proto-cluster gas. 



Table 6. Shock Properties 



T2 


Ti 


Mach No. 


V, 


Compression 


(keV) 


(keV) 


M 


(kms-') 


C 


3.6 ±0.2 


1.5 + 0.4 


2.2 ± 0.5 


1410 ± 110 


2.5 ±0.4 



4. 3. Possible Shock Front Associated with the Radio Relic 

The Coma radio relic, first discovered by Giovannini et al. 
(1985), is one of the best-known rehcs. Its location corre- 
sponds to a projected distance of ~2.1 Mpc from the cluster 
center. The relic measures roughly 800 x 400 kpc in the tangen- 
tial and radial directions, respectively. Thierbach et al. (2003) 
reported a spatially integrated spectral index of a = 1.18 +0.02 
between 151 MHz and 4.75 GHz, consistent with previous 
measurements (Giovannini et al. 1991). 

Our Suzaku data shows a relatively steep temperature drop 
from inside the reUc to outside it towards the southwest. This 
suggests that a shock is located somewhere between the middle 
of the relic and the region just to the southwest;unfortunately, 
the limited spatial resolution of Suzaku does not provide 
stronger constraints on the location of the shock front. Shock 
fronts at the outer edges of radio relics have now been seen in 
a number of other merging clusters (e.g., A3667; Finoguenov 
et al. 2010; Akamatsu et al. 2012). Since the temperature is 
high everywhere to the NE of the reUc, this indicates that the 
shock is propagating towards the SW relative to the pre-shock 
gas. Some of the estimated properties of the shock are listed in 
Table 6. We estimate the Mach number M from the Rankine- 
Hugoniot condition for the temperature jump, 

72 _ 5M^-H4M^-3 

Ti ~ 16M2 ' 
where subscripts 1 and 2 denote pre-shock and post-shock val- 
ues, respectively, and the assumed ratio of specific heats is 
y = 5/3. We estimate the post- and pre-shock temperatures 
based on the temperatures in the annuli from 30'-34' and 43'- 
54', respectively. This gives Al = 2.2 + 0.5, which is similar 
to the values estimated for other shocks associated with radio 
relics (e.g., A3667 Finoguenov et al. 2010; Akamatsu et al. 
2012). With these values, the shock speed is = (1410 ± 1 10) 
km s"'. 

With this value for the Mach number, the shock compres- 
sion is predicted to be C = P2/P1 = 2.5 ± 0.4, where p is the 
gas density. One can also estimate the shock compression 
from the jump in the X-ray surface brightness across the relic. 
The surface brightness was derived from the observed spectra. 
Figure 6 shows the X-ray surface brightness (0.5-10 keV) for 
the same regions used to derived the spectra. There is a drop in 
the surface brightness by a factor of ~5.2 near the outer edge 
of the relic. Projection effects and the low spatial resolution 
of Suzaku are Ukely to cause this to be a slight underestimate 
of the surface brightness jump. For a fixed line of sight path 
through the gas, this jump in surface brightness corresponds 
to a shock compression of C = 2.1, including the difference 
in emissivity for the two regions. Given that this is probably a 
slight underestimate, the surface brightness jump is completely 
consistent with the shock compression of C = 2.5 ± 0.4 pre- 
dicted based on the temperature jump. 

Radio relics generally have rather steep radio spectra, which 
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Fig. 6. Unabsorbed X-ray surface brightness (0.5-10 keV) from the 
X-ray spectroscopy of the annuh to the SW of NGC4839. The red 
crosses are the same amiuh used for the spectra, plus an extra annulus 
just beyond the relic which has too few counts for a temperature de- 
termination. The grey cross combines the two annuli beyond the relic 
to improve the signal to noise. The dashed vertical lines show the ap- 
proximate boundaries of the radio relic. The is a drop in the surface 
brightness near the outer boundary of the radio relic. 

are believed to result from radiative losses by the relativistic 
electrons which produce the synchrotron emission. In a num- 
ber of cases, the radio spectrum is relatively flat at the outer 
edge of the relic at the location of the shock, and steepens 
along the direction of flow away from the shock (e.g., A3667: 
Rottgering et al. 1997, Finoguenov et al. 2010; CIZA2242: van 
Weeren et al. 2010). The width of the relic and the distance 
over which the spectrum steepens are consistent with the pre- 
dicted radiative loss rate and the flow velocity in Abell 3667. 
Assuming the strength of the magnetic field is 1-5 yuG, the ra- 
diative loss time of electrons which radiate at 150 MHz would 
be f,ad » 0. 1 -0.3 Gyr (see Fig. 1 .6 in Sarazin 2002); it is shorter 
for higher frequency radio emission. The velocity of the post- 
shock gas relative to the shock is given by vJC ~ 560 km s"'. 
The width of the relic is about 200-400 kpc, and it takes about 
fad ~ 0.5 Gyr for ICM to travel this distance away from the 
shock. Thus, the radiative loss time is likely to be shorter than 
the advection time for the electrons which emit at all radio fre- 
quencies of at least 150 MHz. 

The magnetic field strengths assumed above can be justified 
from upper limits to the non-thermal, inverse Compton (IC) 
X-ray emission emanating from the same population of 
electrons that generate the radio relic. Relativistic electrons 
that scatter photons of the cosmic microwave background up to 
energies 1-10 keV would have Lorentz factors y~(l-3)xlO^ 
and radiative loss times frad ~ 1 - 3 Gyr, longer than the afore- 
mentioned advection time. Hence the electron spectrum at 
these energies should remain the same as that at acceleration. 



and the expected photon index of the IC X-ray emission is 
r = 1.68. Thus we analyze the XIS spectrum of regions 
around the radio relic by including a power-law compo- 
nent with r = 1.68 in addition to the thermal ICM emission; 
apeci^HB + wabs{apecMV/H + powerlawcxB + apecicu + poweWaitiNonthermai). 
For the spectrum of the 34'-38' annular region along the outer 
edge of the relic (Fig. 3(k)), upper limits on the flux to such 
a component is derived with 7.5 x 10"'''erg/cm^/s/arcmin^ 
(0.3-10 keV). Considering the entire radio-emitting re- 
gion with area S = 800 x 400kpc ~ 380arcmin-, we obtain 
2.9 X 10^'^erg/cmVs (0.3-10 keV). While the observed radio 
flux constrains the combination of the magnetic field strength 
and the density of non-thermal electrons, upper limits to the 
IC flux gives an upper limit to the latter alone, allowing us 
to derive a lower bound on the field strength of B > 0.35 pG. 
This can be compared with the result of Z? > 1 .05pG based 
on XMM-Newton data Feretti & Neumann 2006. It is also in 
Une with values of B ~ 0.5pG obtained from Faraday rotation 
measurements for regions at 2 Mpc distance from the cluster 
center (Bonafede et al. 2010). We note that our flux limits 
are also useful for constraining other potential non-thermal 
emission components (e.g. Sarazin & Kempner 2000; Inoue et 
al. 2005; Vannoni et al. 201 1). 

If the relativistic electrons are accelerated from thermal en- 
ergies by Diffusive Shock Acceleration (DSA) and there is 
no dynamical effect of the relativistic particles or magnetic 
field on the shock, then the energy spectrum of the relativis- 
tic electrons depends only on the shock compression, and is 
given by n(E) dE oc E~p dE, where the power-law index is 
p = {C + 2)/{C-l). For our shock compression, this implies 
3.0 + 0.5. If the magnetic field in the relic is roughly constant, 
such a power-law energy spectrum will lead to synchrotron 
emission with a spectrum 5v oc v"", where a - (p - l)/2 = 
1.00 + 0.25. This should apply to the radio spectrum immedi- 
ately behind the shock. Unfortunately, as far as we are aware, 
the most recent radio spectral index maps for the Coma radio 
relic were given by Giovannini et al. (1991), and they do not 
show any clear systematic spectral variations across the relic. 
Thus, it is most useful to compare to the integrated spectral in- 
dex for the relic. As noted above, we expect radiative losses to 
be quite effective in the Coma relic. If the particles lose energy 
quickly and there is a steady-state between acceleration and 
radiative losses, then the particle spectrum steepens by unity, 
/'loss = p + 1 ~ 2.0 + 0.5, and the integrated radio spectrum is 
steeper by one half, cioss - a + 0.5 1.5 + 0.25. This is steeper 
than the observed integrated radio spectrum, a = 1.18 + 0.02, 
although the errors in the predicted spectral index are large. For 
test-particle DSA, the observed radio spectrum would require 
a shock compression of C = 3.21 + 0.06, which is considerably 
higher than is suggested by the temperature or surface bright- 
ness jump in the shock. Obviously, it would be useful to have 
a better spectral index map for the Coma relic, which might al- 
low us to directly determine the spectrum of the most recently 
accelerated electrons. 

Thus, although our observations confirm the presence of a 
shock associated with the Coma radio relic, and the width of 
the relic is consistent with the flow velocity and radiative loss 
time scale behind the shock, and the radio polarization implies 
that the magnetic field is roughly parallel to the long axis of 
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the relic (Thierbach et al. 2003) as expected for a shock, the 
observed radio spectrum appears to be flatter than would be 
the case for the simplest test-particle DSA theory. Although 
the errors in the shock properties we derived are large, we 
note that similar results have been found for several other 
rehcs, including the NW relic in Abell 3667 (Finoguenov et 
al. 2010; Ogrean et al. 2012; Akamatsu & Kawahara 2011). 
and CIZA 2242 (van Weeren et al. 2010). 

However, compared to the simplest, test-particle theory, the 
physics of DSA under conditions that are realistic and appro- 
priate for clusters is certain to be more complicated, with no 
trivial correspondence between the shock Mach number and 
the spectral index of accelerated particles. There are at least 
two possible ways in which the actual spectrum of accelerated 
particles can be harder than the naive expectation from test par- 
ticle theory, as is observed to be the case here for the Coma 
radio relic. 

The first possibiUty concerns nonUnear DSA in cosmic ray 
modified shocks (see e.g. Malkov & Drury 2001 for a review). 
If the efficiency of injection into the DSA process is suffi- 
cientiy high, the pressure due to the accelerated particles can 
exert a nonlinear, dynamical back reaction onto the structure 
of the shock. The particle spectrum in such circumstances is 
no longer a simple power-law, but attains an overall, concave 
shape. The spectrum at lower energies becomes softer than 
the test-particle case as the velocity jump perceived by these 
particles is reduced by the back reaction, while the spectrum 
at higher energies becomes harder because of the increased 
compressibiUty of the downstream gas and consequently en- 
hanced velocity jump. Nevertheless, this is unlikely to provide 
an explanation for the observed facts here, since 1) shocks with 
Mach numbers as low as ~ 2 are not expected to foster DSA 
with efficiencies high enough for nordinear effects to become 
important, and 2) for the conditions relevant to shocks in clus- 
ters, the harder spectrum is expected to occur at energies much 
higher than for the electrons responsible for the radio rehcs 
(e.g. Kang & Jones 2005), that is, the radio spectrum may ac- 
tually be softer than the test particle case if nonlinear effects 
were effective, and 3) the surface brightness profile shows no 
evidence of extra compression that might characterize nonlin- 
ear shocks, and instead the measured compression agrees rea- 
sonably with (albeit being slightly smaller than) that expected 
from the temperature j ump for a pure gas shock. Another prob- 
lem is that nonlinear shocks are more compressive than pure 
gas shocks, and this accounts for much of the change to the 
spectrum of accelerated particles. However, the measured com- 
pression in the Coma relic shock agrees reasonably with (but is 
slightly smaller than) that expected from the temperature jump, 
implying that there is no extra compression. 

The second and more plausible possibility is DSA initiated 
by injection of pre-existing nonthermal particles, rather than by 
that of thermal particles in the ICM (Kang & Ryu 201 1; Kang 
et al. 2012). Before being processed by the low Mach-number 
shock discussed above, the pre-shock gas may already con- 
tain various types of nonthermal particles, such as those ac- 
celerated by accretion shocks with higher Mach numbers that 
presumably occur at larger radii (Ryu et al. 2003), those ac- 
celerated by turbulence in regions away from shocks (Brunetti 
2011), or those generated continuously as secondary electrons 



and positrons from collisions of accelerated protons with ther- 
mal protons of the ICM (Dennison 1980; Miniati et al. 2001). 
For DSA by low Mach number shocks, the injection of such 
pre-existing nonthermal particles can be much more important 
than that of particles from the thermal distribution for which 
the efficiency is expected to be quite low. In such cases, the 
spectrum of accelerated particles may be determined by that of 
the pre-existing population rather than the shock Mach num- 
ber (Kang & Ryu 2011; Kang et al. 2012). Kang et al. 2012 
have specifically shown that such a model can reproduce the 
observed properties of some known radio rehcs associated with 
shocks detected in X-rays. Although a quantitative discussion 
is beyond the scope of this paper, we note that for the Coma ra- 
dio rehc, the pre-existing nonthermal electrons must have been 
produced not far from the shock front to avoid excessive ra- 
diative losses. If they were secondaries resulting from proton- 
proton collisions, the parent protons may have originated fur- 
ther in the past, for example from within the NGC 4839 group 
before its passage through the virial radius, or may have been 
advected inward along with cooler gas from the large-scale 
structure (LSS) filament in the SW direction, or both 

4.4. Origin of the Shock Front at the Radio Relic 

In modern theories of hierarchical structure formation, both 
internal merger shocks and extemal accretion shocks around 
clusters of galaxies are predicted (Miniati et al. 2000; Ryu et 
al. 2003). External accretion shocks are located at large radii 
(about 3 times the virial radius), are more spherical in shape, 
and are expected to have high Mach numbers of up to 10^. 
Internal merger shocks located around the virial radius have 
Mach numbers ranging between 2-10 in simulations, and are 
smaller in extent. The temperature jump associated with the 
Coma radio relic located at 2 Mpc (0.9 r2oo) gives a Mach 
number of Al ~ 2, which is not expected for an extemal accre- 
tion shock from the diffuse intergalactic medium. This might 
be an outgoing merger shock associated with the infall of the 
NGC4839 group. However, the tail and abundance structure 
around NGC4839 suggest that it is moving to the NE, so the 
relic shock is not a merger bow shock in front of this group. 
Another possibility is that the rehc shock is due to accre- 
tion from the large scale structure filament (Brown & Rudnick 
2011) which extends from Coma SW towards the Abell 1367 
cluster. Recently, Kale et al. (2012) discuss the effect of large- 
scale accretion flows on the outgoing merger shocks based on 
their radio observations of the typical major merger cluster 
Abell 3376. Recent numerical simulations (Paul et al. 2011) 
indicate that shock fronts are enhanced along filament direc- 
tions. In addition, the radio morphology of the radio galaxy 
NGC4789 suggests that it is affected by an accretion flow onto 
the Coma cluster from the large scale structure filament (see 
Fig. 1 in Giovannini et al. 1991). Considering the lower ICM 
temperature at the pre-shock region and the evidence for the 
infall of NGC4839, the shock is may be generated by the com- 
bination of an outgoing merger shock and the inflow of cooler 
gas from the large scale structure filament towards Abell 1367. 

Our data may provide direct observational evidence of a 
shock involved with the infall of the NGC4839 group. Such 
shocks are a means to transform the kinetic energy of the merg- 
ing subcluster into the thermal energy of the ICM, and this pro- 
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cess can affect the thermal history of clusters of galaxies. As 
shown in our previous work (Akamatsu & Kawahara 20 11), ra- 
dio rehcs are good tracers of large scale structure shock fronts. 
Further studies (radio emission, SZ, X-rays, and simulations) 
are needed in order to reveal the nature of these shock fronts . 

5. Summary 

We observed NGC4839 and the Coma radio relic with the 
Suzaku XIS and derived radial profiles of temperature and 
abundance. The temperature across radio reUc shows a jump, 
which indicates the presence of a shock front. We estimated a 
Mach number of Al ~ 2.2 using the Rankine-Hugoniot temper- 
ature jump condition. 

The main results of this work are sunomarized as follows; 

• The ICM temperature distribution is flat between 10' - 

30' from NGC4839, and then decreases sharply (consis- 
tent with a M ~ 2.2 shock) at some location near the 
outer edge of the radio rehc. 

• Compared to the inner regions, there is a clear change in 
the temperature structure around NGC4839. 

• The ICM abundance has a peak at NGC4839 and grad- 
ually decreases from 0.6 to 0.2 solar towards the outer 
regions. The abundance profile shows an excess over the 
mean values for other clusters at such large radii as ob- 
served with XMM-Newton. 

• The estimated Mach number of the shock is Al=2.2±0.5 
based on the temperature jump, which corresponds to 
a shock velocity of v, = (1410+ 110) km s^'. The X- 
ray surface brightness also shows a drop near the outer 
edge of the radio relic. Applying the shock density jump 
condition to the surface brightness profile gives a Mach 
number which is consistent with the value from the tem- 
perature jump. The implied post-shock velocity and the 
rate of radiative losses by the radio-emitting relativistic 
electrons imply that losses are very significant, and can 
help to explain the relatively narrow width of the relic. 

• The observed integrated radio spectrum of the reUc is 
somewhat flatter than would be predicted by test-particle 
diffusive shock acceleration (DSA) in a shock with the 
observed Mach number and compression. A similar re- 
sult is seen in other relics. We discuss possible expla- 
nations for this discrepancy, and conclude that the shock 
may be re-accelerating an existing population of low en- 
ergy relativistic electrons, rather than directly accelerat- 
ing thermal electrons to high energies. This might also 
explain why some shocks are seen in clusters without as- 
sociated radio relics. 

These results show that the NGC 4839 group is falling into the 
Coma cluster, and that the radio relic was possibly generated 
by shock acceleration. The origin of the shock structure is still 
unclear Upcoming radio (LOFAR: van Weeren et al. 2012) and 
SZ (ALMA: Yamada et al. 2012) observations will provide new 
insights about the nature of the shock front and the dynamical 
history of the Coma cluster. Obviously, additional wide field 
X-ray observations would also be very useful. 

As this work was being completed, we became aware of the 
preprint by Ogrean & Briiggen (2012), who report evidence for 



a shock with Mach number Al ~ 2 at the Coma radio relic based 
on archival XMM-Newton data. 
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